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Abstract 

production in two-photon collisions has been studied using a large 
data sample of 67 fb~' accumulated with the Belle detector at the KEKB 
asymmetric e"'"e~ collider. We have measured the cross section for the pro- 
cess 77 K~^K~ for center-of-mass energies between 1.4 and 2.4 GeV, and 
found three new resonant structures in the energy region between 1.6 and 
2.4 GeV. The angular differential cross sections have also been measured. 



PACS. 13.20.Gd Two-photon collisions - 13.60.Le Meson production 
- 14.40.Cs Other mesons with S = C = 
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1 Introduction 



A high luminosity electron-positron collider is well suited for studies of meson reso- 
nances produced in two-photon collisions. The heaviest established resonance that 
has been observed so far in kaon-pair production in two-photon processes is the 
/2(1525) meson, which is classified as an almost pure ss meson. Above the /2(1525) 
mass, no clear resonance structure has been found in the K^K^ channel [1, 2]. 

The L3 experiment at LEP has reported a resonance- like peak at 1.767 GeV in 
the 77 — > K^K^ process [3]. Their recent analysis shows a dominant contribution 
of a tensor component in this energy region. However, assignment of this structure 
to any known resonance state or a spin/isospin state is not yet conclusive. The 
corresponding resonant structure is expected to appear also in the K^K~ channel 
because of the isospin invariancc of resonance decay. However, we cannot predict 
its shape in the cross section since two-photon reactions themselves are not isospin 
invariant and the continuum contributions and interference effects can differ between 
the K^K^ and K'^K~ channels. Therefore, measurements of the 77 K'^K~ 
and 77 K'gK^ processes in the same mass region give essentially independent 
information. 

The 1.6 - 2.4 GeV region is very important for meson spectroscopy since radially- 
excited qq states are expected to exist in this mass range. Several mesons with 
poorly measured properties are reported in this region having J^^ = (even)"'"'' [4], 
where J, P and C are spin, parity and charge conjugation, respectively. The prop- 
erties of these mesons can be measured precisely using a large clean sample of 
77 KK events. 

In addition, two- photon processes play an important role in identifying glueballs. 
since we expect much smaller coupling of photons with a glueball than with a qq 
meson. Some glueball candidates around 2 GeV have been observed inJ/tp^ -jKK 
decays [5], but none has been identified conclusively as a glueball. 

Here, we report results on measurements of 77 K~^K~ for two-photon center- 
of-mass energy between 1.4 and 2.4 GeV. We describe the experimental apparatus 
and triggers in section 2. In section 3, the selection criteria of the signal events are 
shown. Section 4 introduces the methods used to derive the total and differential 
cross sections and presents the results. The major sources of systematic errors in the 
present measurement are itemized in section 5. We apply phenomcnological analyses 
to the obtained cross sections, mainly with respect to extractions of resonances 
in this energy region, in section 6, and discuss the attributes of these resonant 
structures in section 7. Conclusions are drawn in section 8. 

2 Experimental Data and the Belle detector 

The experiment is carried out with the Belle detector [6] at the KEKB asymmet- 
ric e+e" collider [7]. In KEKB, the electron beam (8 GeV) and positron beam 
(3.5 GeV) collide with a crossing angle of 22 mrad. The data collected between 
1999 and April 2002 are used for this analysis, corresponding to an integrated lu- 
minosity of 67 fb~^. We combine data taken at the on- and off-resonance energies; 
the off-resonance data are taken 60 MeV below the T(45) resonance at 10.58 GeV. 

In the present analysis, neither the recoil electron nor positron is detected. The 
basic topology of the signal events is just two charged tracks with opposite charge. 
We restrict the virtuality of the incident photons to be small by imposing a strict 
transverse-momentum balance rcqiiirement on this two-track system with respect 
to the incident axis in the e^e^ center of mass (cm.) frame. 

A comprehensive description of the Belle detector is given elsewhere [6]. We 
mention here only the detector components essential for the present measurement. 
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Charged tracks are reconstructed from hit information in a central drift chamber 
(CDC) located in a uniform 1.5 T solenoidal magnetic field. The z axis of the 
detector and the solenoid are along the positron beam, with the positrons moving 
in the —z direction. The CDC measures the longitudinal and transverse momentum 
components (along the z axis and in the r</3 plane, respectively). The transverse 
momentum resolution is determined from cosmic rays and e+e~ events 
to be {(Jpi/ptY = 0.0030^ + (O.OOlQpt)^, where pt is the transverse momentum 
in GeV/c. Track trajectory coordinates near the collision point are provided by 
a silicon vertex detector (SVD). Photon detection and energy measurements are 
performed with a Csl electromagnetic calorimeter (ECL). In this analysis, the ECL 
is mainly used for rejection of electrons. Identification of kaons is made using 
information from the time-of-flight counters (TOP) and silica-aerogel Cherenkov 
counters (ACC). The ACC provides good separation between kaons and pions or 
muons at momenta above 1.2 GeV/c. The TOP system consists of a barrel of 128 
plastic scintillation counters, and is effective in K/ir separation mainly for tracks 
with momentum below 1.2 GeV/c. The lower energy kaons arc identified also using 
specific ionization {dE/dx) measurements in the CDC. The magnet return yoke 
is instrumented to form the and muon detector (KLM), which detects muon 
tracks and provides trigger signals. In this analysis, the KLM is used only for the 
calibration of the trigger efficiencies using muon-pair events; it is not needed nor 
used for particle identification of the signal events. 

Signal events are triggered most effectively by requiring two or more CDC tracks 
in the r^p plane as well as two or more TOP hits and at least one isolated cluster 
in the ECL with an energy above 0.1 GeV. The opening angle in the r^p plane of 
the two tracks must exceed 135°, and at least one of them must have z coordinate 
information from the CDC cathodes. Additional triggers that combine the two-or- 
more track requirement with either a KLM track or an ECL energy deposit above 
0.5 GeV also collect some signal events; they are used in conjunction with other 
redundant triggers to detect muon-pair and electron-pair events from two-photon 
collisions and thereby calibrate the trigger efficiency for our signal events. 

The trigger efficiency for the K^K~ process is 80-93% in most of our selection 
acceptance, with the variation depending on the transverse momentum of the tracks. 

3 Event Selection 

Candidate events arc selected offline using the following criteria. The event must 
have two oppositely charged tracks, each with pt > 0.4 GeV/c, \dr\ < 1 cm, \dz\ < 
3 cm, and —0.47 < cos6' < -1-0.82, and both satisfying \dzi — dz2\ < 1 cm. Here, 
Pt is the transverse momentum with respect to the positron-beam axis, \dr\ and dz 
are the radial and axial coordinates, respectively, of the point of closest approach 
to the nominal collision point (as seen in the rip plane), and 9 is the polar angle, all 
measured in the laboratory frame. We demand that the event have no extra charged 
track with pt > 0.1 GeV/c, and then use only the two tracks with pt > 0.4 GeV/c 
in this analysis. The sum of the track momentum magnitudes must be less than 
6 GeV/c, and the total energy deposition in the ECL must be less than 6 GeV. 

Events with an initial-state radiation, such as radiative Bhabha events, are sup- 
pressed by requiring that the invariant mass of the two tracks be smaller than 
4.5 GeV/c^ and the square of the missing mass of the event be greater than 
2 GeV^/c^. Cosmic ray events are rejected by requiring the cosine of the open- 
ing angle of the tracks to be greater than —0.997. Exclusive two-track events from 
quasi-real two-photon collisions are selected by requiring a good transverse mo- 
mentum balance in the e+e~ cm. frame for the two energetic tracks: | X^Pj | = 
|Pn+Pt*2l <0.1 GeV/c. 
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After the application of these selection criteria, 2.56 x lO'^ events rciniain. They 
are dominated by two-photon events of light-particle pairs: 77 — > e+e", /x+/Lt~ and 

TT+TT". 

We apply particle identification to the tracks in the remaining events to select 
77 K^K~ events. We suppress e+e~ events by requiring E/p < 0.8 for each 
track (where E is the energy of the ECL cluster that matches the track of momentum 
p). We require that the TOF counter system gives useful time-of-flight information 
for each track. Wc suppress TT^n~, jJ,'^ ix~ and pp events by requiring that the 
TOF-, ACC- and (for K/p separation only) dE / dx-dev\ved particle identification 
likelihood ratios Ck/{Ck + C-j^) and Ck / {Ck + Cp) each exceed 0.8 for each track, 
and demanding ATOFk > — 0.6 ns for each track, where ATOFk is the difference 
between the TOF-measured and calculated times of flight of the charged track, 
assuming it to be a kaon. 

An ambiguity in the determination of the collision time {to) introduces spurious 
events into our sample at low invariant mass {Mj^+k- < 1-55 GeV/c^). Our to 
measurement in each event initially has an ambiguity among times that differ by 
an integer multiple of the RF-beam-bucket spacing time {tg = 1.965 ns), and we 
determine this off'set in each event by requiring the consistent TOF identification 
of both tracks as a pair of known long-lived charged particles. However, there is a 
kinematical region where the to assignment has a two- fold solution, differing by ±ts, 
that can lead to misidentification of 7r"'"7r~//i+/i^/e+e~ as K~^K~ events or vice- 
versa. To suppress the contamination from the lighter partic;le-pair events where 
the wrong to was chosen, we apply the following prescription. Assuming that the 
two energetic tracks are pions, we calculate ATOF^ as well as xi = (^T'OF^/dt)^, 
where at = 0.10 ns is the nominal time resolution of the TOF, for each track using a 
collision time that is either nominal or shifted forward or backward by tg, then select 
the case with the minimum + X^- • We do the same, assuming both tracks to 
be kaons. If the optimal collision time for the kaon hypothesis is one beam bucket 
interval earlier than the optimal time for the pion hypothesis, then we reject the 
event if xl+ + xl- < 49 or ATOF^+ + ATOF^- > -0.73 ns. If the event survives 
this test, we also require x'k+ + Xk- < ^5 and that the dE/dx measurements in 
the CDC for each track be consistent with a kaon. 

Proton backgrounds and the loss of signal due to KK misidentified as pp because 
of such to misassignment are confirmed to be negligibly small by using a tighter 
dE/dx cut and by investigating the ATOFk distribution, respectively. 

Figure 1 shows the ATOFk distribution of a track in an event where the other 
track is positively identified as a kaon and where both tracks are so identified. 
(The latter samples correspond to the final signal candidates.) The true kaons — 
the peak near ATOFk = ns — are well separated from lighter particles (the peak 
below zero). Clearly, the lighter-particle backgrounds are rejected while nearly 
all K^K~ events are retained by these cuts. The largest remaining background 
contribution from particle misidentification comes from 77 — > /2(1270) — > tt+tt", 
which contaminates the signal sample up to 6% in a narrow K~^K~ invariant mass 
region around 1.60 GeV. 

After the application of all the selection criteria, 63455 candidate events remain. 
Their invariant-mass distribution is shown in Fig. 2. A peak at aroimd 1.52 GeV 
comes from the /2(1525) resonance. A decrease of the events below 1.45 GeV is 
mainly due to an effect from the selection criteria for the pt cut of tracks and the 
cut to avoid the to determination ambiguity. A few bump structures are seen at 
higher energies. 
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Figure 1: The ATOFk distribution for each track in the event where the other 
track is identified as a kaon (histogram), and where both tracks are so identified 
(asterisks), in tliree different W regions. 
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Figure 2: The invariant mass distribution of the K^K candidates that satisfy all 
the selection criteria. 
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4 Derivation of Cross Sections 



The cross section for 77 — > K~^K~ is derived from the present measurement. We 

restrict the polar-angle range of the final state kaons in the 77 cm. frame (9*) to 
be within | cos9*\ < 0.6. The differential cross sections are given by the following 
formula: 

d<T ^ AN{W,\cose*\) (1-/bgTO ^ 

d\cos0*\ jCdt[Lj^{W)AW] A\cose*\r){W,\cose*\)' ^' 

where AN{W, | cos^*|) is the number of events in a two-dimensional bin of the cm. 

energy (W) and the cosine of the polar angle of a meson in the cm. system (| cos 0*\) 
with their widths AW and A|cos0*|, respectively, where we take the absolute 
value of the cosine according to the symmetry of the initial-state particles. In 
the denominator, / Cdt, _L^.-j,(Ty), and rj(W, \ cos6'*|) arc the integrated luminosity, 
the luminosity function and the event efficiency at the given kinematical point, 
respectively. 

Residual background due to particle misidcntification is corrected using an 
energy-dependent factor fBG{W), that is derived from a study of the ATOFk 
distributions shown in Fig. 1. This background fraction is estimated to be less than 
3% except in the range 1.55 - 1.65 GeV where it can be as large as 6% due to 
/2(1270) TT+TT^ contamination. We neglect any potential 9* dependence of the 
correction factor since no prominent angular variation is anticipated in the back- 
ground that arises mainly from the non-Gaussian tails of the TOF measurements. 

The estimated contribution from other background processes, discussed in sec- 
tion 5.3, is much smaller than the total systematic errors (shown in Table 2). We 
neglect the effect of such backgrounds in the derivation of the cross section. 

We use the measured invariant mass of the K'^K~ system as the 77 cm. energy 
W. The fractional energy resolution AW/W ci: 0.2% is estimated from a Monte 
Carlo (MC) simulation of the signal process (described below). Since this is much 
smaller than our energy bin size, we neglect smearing across bins in the cross section 
derivation. 

Since we cannot determine the true 77 axis in each event, we instead measure 
9* from the direction of the e^e'^ beam axis in the e^e~ cm. frame. The difference 
between this and the true polar angle is confirmed to be small using the signal MC 
simulation, corresponding to an r.m.s. deviation of about 0.015 in cos9* . 

The efficiency is determined for each bin using the signal MC events for the 
e+e~ e~^e~K~^K~ process that are generated by TREPS [8] and simulated 
within the Belle detector by a program based on GEANT3 [9]. This efficiency is 
represented by a smooth function of {W, \ cos 6**1) in the analysis. 

The trigger efficiency for K~^K~ events is evaluated separately via a study of 
detected two-photon e~^e~ and events that satisfy two or more independent 

trigger conditions. We parameterize the trigger efficiency as a function of the av- 
eraged transverse momentum (pt) of the two tracks in an event. The transverse 
momentum difference is at most 0.1 GeV/c because of the selection criteria; under 
these circumstances, the relation pt — \J {W/2y — tuk^ sin9* is quite accurate, so 
the trigger efficiency is calculated by this formula in each {W , \ cos9*\) bin. The 
laboratory-angle dependence of the trigger efficiency is small within the acceptance 
range, and we neglect the effect. 

The typical value of the efficiency rj{W, \ cos 9* \ ) is 13% at large angles ( | cos 9* \ < 
0.3) in the high-invariant mass region (W > 1.8 GeV). At most of the other points, 
it is between 1% and 10%. The large inefficiency seen in rj(W, | cos 6**1) is mainly 
attributed to detector acceptance effects. The efficiency to identify kaon-pair events 
within the acceptance is typically 55%, including losses from kaon decays. 



7 



We calculate the luminosity function using a separate feature of TREPS [8]. 
The effects of longitudinal photons are ignored therein. We introduce the form 
factor effect for finite-Q^ photons by multiplying the two-photon flux by the factor 
{1 + Ql/W'^)~'^{l + Q2/W^)~'^ in the integrations over and Q2 for the calculation 
of the luminosity functions. (Here, and Q2 are the absolute values of the four- 
momentum transfer of each photon.) The tight cut of 0.1 GeV/c on | Pt I this 
analysis limits the uncertainty from the choice of the form factor to less than 2%. 
The systematic error of the luminosity function is estimated to be 4% by a study that 
compares the cross sections for the e+e" — > e~^e~ fi~^ iJ,~ process using this luminosity 
function to those from a full a^-ordcr QED calculation [10]. The luminosity function 
L^j{W) and the efficiency r]{W, \ cos 6**1) are consistently calculated for two-photon 
collisions in the same Ql and Q2 ranges, Ql, Q\ <1 GeV^. The calculated Lry^{W) 
is a smoothly decreasing function equal to 5.9 x 10~^ GeV""*^ aiW = 1.4 GeV and 

2.0 X 10-3 GeV-i aiW = 2.4 GeV. 

The cross section values integrated over | cos^*| < 0.6 are derived by a binned 
maximum likelihood fit of the experimental | cos 9* \ distribution with a bin width of 
0.05 carried out at each W bin with a bin width of 0.02 GeV. We assume that the 
angular dependence of the differential cross section is a second-order polynomial of 
cos^ 9* and integrate the fit with | cos6'*| between and 0.6. Several data points at 
W < 1.46 GeV and | cos9*\ > 0.50 are removed from the fit, where the efficiencies 
are extremely small. 

The obtained cross sections for the process 77 K^K~ integrated over the 
range | cos^*| < 0.6 are summarized in Table 1 and depicted in Fig. 3. The error 
bars in the figure are statistical only. Systematic errors are tabulated in Table 2, 
and the dominant sources are described in section 5. 

Figure 3 shows our result in comparison with the previous measurements. The 
errors for our data are statistical only; we note that they are much improved over 
those of ARGUS [1] and TPC/Two-Gamma [2]. The ARGUS cross sections are 
obtained for the full angular range (|cos^*| < 1) by fitting the angular distribu- 
tions to the sum of two spin-helicity components, (J, A) = (0,0) and (2,2), and 
should therefore be somewhat larger typically by 30% — than those of our mea- 
surement and TPC/Two-Gamma. (See section 6 for the spin-helicity definition and 
decomposition.) 

Our differential cross sections in the 77 cm. polar angle arc plotted in Fig. 4 for 
each energy bin of width AW = 0.04 GeV. Again, the displayed errors are statistical 
only. 

5 Major sources of systematic error 

5.1 Trigger efficiency 

We determine the trigger efficiency for the final samples that would pass all 

the selection criteria described in section 3. We first estimate the efficiency of the 
two-track trigger, which requires two or more TOF/TSC hits as well as one or more 
Csl clusters. The efficiency is determined experimentally as a function of transverse 
momentum between 0.4 and 1.2 GeV/c, using the redundancy of triggers for two- 
photon e^e~ and events. The ratio of the yield from this particular trigger to 
that from all the trigger sources — typically 0.88 — is used to determine the overall 
trigger efficiency. We confirm that there is no notable explicit dependence of the 
trigger efficiency on the polar angle in the laboratory system. The results on the 
trigger efficiency are compared with those from our trigger simulation program, 
which is applied immediately after the full detector simulation. The simulation 
reproduces a quantitative nature of the pt dependence of the trigger efficiency. 
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Table 1: The cross section of the process 77 K~^K in the polar angular region 
I cos^*| < 0.6. The first error is statistical and the second systematic. 



l^(GeV) 0- (|cosr| < 0.6) (nb) 



1, 


.40 - 


1 


.42 


5.7 ±0.7 ±1.2 


1, 


.42 - 


1 


.44 


5.4 ± 0.7 ± 1.0 


1, 


.44 - 


1 


.46 


5.4 ±0.5 ±0.8 


1 


.46 - 


1 


.48 


7.1 ± 0.3 ± 1.0 


1, 


.48 - 


1 


.50 


10.0 ± 0.2 ± 1.4 


1, 


.50 - 


1 


.52 


14.4 ±0.3 ±2.0 


1 


.52 - 


1, 


.54 


14.2 ± 0.3 ± 1.9 


1, 


.54 - 


1 


.56 


11.4 ± 0.2 ± 1.6 


1, 


.56 - 


1 


.58 


7.9 ± 0.2 ± 1.0 


1, 


.58 - 


1, 


.60 


5.98 ±0.15 ±0.77 


1, 


.60 - 


1 


.62 


4.81 ±0.13 ±0.64 


1, 


.62 - 


1 


.64 


4.22 ±0.12 ±0.56 


f. 


.64 - 


1 


.66 


3.70 ±0.11 ±0.48 


1, 


.66 - 


1 


.68 


3.55 ±0.11 ±0.44 


1, 


.68 - 


1 


.70 


3.31 ±0.11 ±0.41 


1 


.70 - 


1, 


.72 


3.27 ±0.10 ±0.41 


1, 


.72 - 


1 


.74 


3.29 ±0.11 ±0.41 


1, 


.74 - 


1 


.76 


3.23 ±0.10 ±0.40 


f. 


.76 - 


1 


.78 


3.00 ±0.10 ±0.37 


1, 


.78 - 


1 


.80 


2.94 ±0.10 ±0.36 


1, 


.80 - 


1 


.82 


2.62 ±0.09 ±0.32 


1 


.82 - 


1 


.84 


2.71 ±0.09 ±0.32 


1, 


.84 - 


1 


.86 


2.56 ±0.09 ±0.30 


1, 


.86 - 


1 


.88 


2.64 ±0.09 ±0.31 


1 


.88 - 


1 


.90 


2.92 ±0.09 ±0.34 


1, 


.90 - 


1 


.92 


2.95 ±0.09 ±0.35 


1, 


.92 - 


1 


.94 


2.63 ±0.08 ±0.31 


1, 


.94 - 


1 


.96 


2.84 ±0.09 ±0.33 


1, 


.96 - 


1 


.98 


2.92 ±0.10 ±0.34 


1, 


.98 - 


2.00 


2.59 ±0.09 ±0.31 


2, 


.00 - 


2.02 


2.26 ±0.08 ±0.26 


2, 


.02 - 


2.04 


2.36 ±0.08 ±0.27 


2 


.04 - 


2.06 


1.96 ±0.08 ±0.23 


2, 


.06 - 


2 


.08 


1.68 ±0.07 ±0.20 


2, 


.08 - 


2 


.10 


1.57 ±0.07 ±0.18 


2, 


.10 - 


2 


.12 


1.34 ±0.06 ±0.17 


2, 


.12 - 


2, 


.14 


1.15 ± 0.06 ± 0.14 


2, 


.14 - 


2 


.16 


0.96 ±0.06 ±0.12 


2, 


.16 - 


2 


.18 


0.89 ±0.06 ±0.11 


2, 


.18 - 


2 


.20 


0.89 ±0.06 ±0.11 


2, 


.20 - 


2 


.22 


0.97 ±0.06 ±0.12 


2 


.22 - 


2 


.24 


0.90 ±0.06 ±0.11 


2, 


.24 - 


2 


.26 


0.83 ±0.06 ±0.10 


2, 


.26 - 


2, 


.28 


0.84 ±0.06 ±0.10 


2, 


.28 - 


2.30 


0.96 ±0.06 ±0.12 


2, 


.30 - 


2.32 


0.80 ±0.05 ±0.10 


2, 


.32 - 


2.34 


0.93 ±0.06 ±0.12 


2, 


.34 - 


2.36 


0.76 ±0.06 ±0.09 


2, 


.36 - 


2.38 


0.79 ±0.06 ±0.10 


2, 


.38 - 


2.40 


0.68 ±0.06 ±0.08 
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Table 2: The systematic errors for the cross section 77 K~^K~ in the polar 
angular region | cos^*| < 0.6. Some systematic errors are H^-dependent; these are 
shown as ranges. 



Source 


Systcmiatic 
error 


Trigger efficiency 


4 


-6% 


Tracking efficiency 




4% 


E/p cut 


4 


-6% 


TOF efficiency 




3% 


Kaon identification efficiency by ACC 





- 5% 


to ambiguity 





- 7% 


Acceptance calculation 


3 


- 15% 


Particle misidentification backgrounds 


1 


- 2% 


Non-exclusive backgrounds 


4 


-6% 


Integrated luminosity 




1% 


Liuninosity fiurction and form-factor (effect 




4% 


Total 


11 


- 20% 
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Figure 3: The cross section for 77 — > K~^K~ in the polar angular region 
I cos 6** I < 0.6 obtained in this experiment compared with results from previous 
experiments [1,2]. The error bars are statistical only. 
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Figure 4: The differential cross section da/d\ cos9*\ for 77 K^K~ obtained in 
this experiment. The sohd curves are the fits described in section 6.2. The error 
bars are statistical only. 
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Although the trigger cfRcicncy value from the simulation is not used for the cross 
section calculations, the difference between experimental and MC values is taken 
into account in the systematic error. A small difference between the trigger efficiency 
for leptons and hadrons estimated from the experimental and simulation studies is 
applied as a correction and included to the systematic error as well. The trigger 
efficiency thus determined is (82 ± 5)% at pt = 0.55 GeV/c and (92 ± 4)% at 
Pt = 1.0 GeV/c. 

The samples after application of the event selection criteria before particle iden- 
tification are used for a confirmation of the trigger efficiency. We have compared the 
experimental yields of these samples in the W and laboratory polar angle distribu- 
tions to the expectations from the MC calculation. Those samples are dominated 
by e+e~ e+e~e+e~, e+e~^"'"yu~ and e+e~7r+7r~ processes whose cross sections 
are well known [10, 11, 12]. The experimental and MC yields are consistent within 
the above error of the trigger efficiency. 

5.2 Kaon-identification efficiency 

We have checked the kaon-identification efficiency in this analysis based on the real 
data. The efficiency to identify and select kaon-pair events within the detector 
acceptance is typically 55%. The 45% loss is partially due to kaon decays, which 
account for 20-30% depending on the two kaon invariant mass. An additional 
reduction of ^ 25% comes from the kaon identification criteria for the TOF, ACC 
and E/p requirements. 

The efficiency of the TOF counters giving a useful TOF measurement for a 
track, which is required in the kaon selection, is obtained using the events in which 
the other track is tagged as a kaon. The efficiency thus obtained in the real data is 
around 92% in the 0.5 < p < 1.0 GeV/c region where the signal events dominate, 
comparable to the MC simulation's efficiency of 94% at around 1 GeV/c. The effect 
of a trigger efficiency loss due to kaon decays is small in this comparison.^ We take 
the difference of the TOF efficiency between the experimental and MC data, 3% for 
the two tracks, as a systematic error. 

The efficiency of kaon-pair identification in the offiine selection is examined using 
the redundancy of two kaons in a signal K^K~ event. The ATOFk distribution 
for a track is investigated for events where the other track is identified as a kaon 
as shown in Fig. 1. We find that charged tracks whose ATOFk is consistent with 
a kaon are identified as a kaon with an efficiency higher than 95%; this value is 
consistent with the MC expectation within the systematic error of 5% or less in all 
W regions. 

The E/p cut gives an inefficiency in kaon selection as a nuclear interaction of a 
negative kaon in the ECL sometimes gives a large energy deposition. We compared 
the fractions of kaons with E/p> 0.8 to those with any E/p between the MC events 
and real K~^K~ samples; they are consistent. A small possible discrepancy (~ 6% 
a.t W < 1.8 GeV and ~ 4% at W > 1.8 GeV) is accounted for in the systematic 
error. 

5.3 Backgrounds from pEirticle misidentification and non- 
exclusive events 

The correction factor for the background contamination from particle misidentifica- 
tion is evaluated from the study of the experimental data, as described in section 4. 

'^The effect of kaon decays in flight is taken into account in the simulation, and it appears as 
a loss of the TOF efficiency (and thus, the loss of the K~ signals) in both experimental and 
simulated data. However, a few of the decays in the experimental data also induce a loss at the 
trigger stage; direct comparison is impossible in the low-energy region. 
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Figure 5: The transverse momentum balance (|X]Ptl) distribution for the signal 
process from the real data (closed circles) and signal MC events (histogram) for 
events in three W regions. The MC distributions are scaled to match the data yield 
in the first two bins. 



Another study using the MC expectations from the e~^e~e'^e~ , e~^e^ fi^fi^ [10] and 
e+e^TT+TT^ processes [11, 12], which dominate this kind of background, is also con- 
sistent with this evaluation. We assign a systematic error of 1-2% from this source, 
depending on W. This background source has no effect on the significance of the 
resonant structures that are observed in the spectrum and discussed in this 

paper. 

We also estimate the background contamination from events where additional 
particles accompany the two detected tracks — so-called non-exclusive backgrounds. 
Such events are expected to give a larger | X] Pt I than the signal events. Figure 5 
shows the distribution of | X] Pt I the K'^K~ signal candidates, with the sim- 
ulated distribution scaled to match the data yield in the first two bins. The data 
and signal MC show generally good agreement over the entire accepted range of 
I X Pt 1 7 s-nd we conclude that the backgrounds from non-exclusive processes are 
small. The small difference between the experimental and MC distributions seen 
at high I X Pt I '^^y tie attributed to the non-exclusive backgrounds (or any other 
process that does not have an enhancement at zero-pt balance). 

We find that the contamination of such background processes is less than 6% 
at any W , based on the expectation that the background contribution vanishes at 
I X Pt I = and the assumption that the observed surplus over MC expectation in 
the rightmost bin of Fig. 5 is purely from the contribution from the non-exclusive 
background. In the experimental data, a larger excess over the MC expectation 
is seen at jXPtl around 0.01 - 0.04 GeV/c. This excess cannot be attributed 
to the non-exclusive backgrounds, and is considered to be due to an unmodelled 
broadening of the signal-process distribution.^ 

It is difficult to determine the magnitude of these backgrounds, so we neglect 

^This broadening may be due to interferences or higher-order processes in the relevant photon- 
emission diagrams, as well as a mismatch in the real and simulated resolution of the detector 
components. 
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their effect in the derivation of cross sections and instead account for them in the 
systematic error. 



6 Phenomenological analyses 
6.1 The c.m.-energy dependence 

The obtained 77 c.m.-energy dependence of the experimental yield (Fig. 2) and 
the cross section (Fig. 3) clearly show four resonant structures, in the vicinities of 
1.52 GeV, 1.7 - 1.8 GeV, 1.9 - 2.1 GeV, and 2.2 - 2.4 GeV. The peak near 1.52 GeV is 
from the /2(1525) resonance, which is considered to be a groimd-statc tensor meson 
dominated by an ss component. The contribution of this resonance in this reaction 
is well known from previous measurements [1, 2]. In contrast, this measurement 
constitutes the first observation of the other three structures. 

The W dependence of the cross section is fitted to a sum of contributions from 
resonances and a continuum component. The components adopted in the fit are: (A) 
the /2(1525) resonance, (B) a resonance around 1.75 GeV, (C) a resonance around 
2.0 GeV, (D) a resonance around 2.3 GeV, (E) contributions from the continuum 
component and the tails of two low-mass resonances, /2(1270) and 02 (1320). For 
(A)-(D), we use an amplitude with a relativistic Breit-Wigner form: 

" T4^2_M2-iMr' 

where 5, M and F are the size parameter, mass and total width of the resonance, 
respectively. The component (E) is expressed by an amplitude with a W-dependent 
form and three free parameters a, b, and Wq as 

where we combine the effects from the tails of the low-mass resonances and the 

smooth continuum since we cannot decompose these effects in our present mea- 
surement. The interference phases between (A) and (E) and between (C) and 
(E) are treated as free parameters. (B) and (D) are considered incoherent since 
their peaks are relatively small and it is difficult to determine the rela- 
tive phase. For the /2(1525) resonance, the W-dependent total width is 
T{W) = Vm{p* IPlf{M/W){D2{p*r)/D2{plr)), where M and Tm are the nomi- 
nal mass and total width, respectively, while p* and Pq are the momenta of the 
final-state kaon in the rest frame of the resonance with invariant mass W and M, 
respectively. The function D2{z) ~ (9 -|- iz^ -\- z'^)~^ [13] represents a centrifugal 
barrier factor; we use an effective interaction radius of r = 1 fm. We use W- 
indcpendent total widths for the other three resonances. We denote the Fm for the 
/2(1525) simply by F in the remainder of this paper. 

The result of the fit is summarized in Table 3. The curves in Fig. 6 show the 
fit and the contribution from each component superimposed on the experimental 
data. The goodness of fit is x^/ndf = 50.8/33. We use only statistical errors for 
this calculation and its minimization. We determine the systematic errors of the 
fitted parameters by a study of their changes for four cases of linear deformations 
of the measured cross section as a function of 



^Moreover, it is noted that the interference term from two distinet resonance amplitudes van- 
ishes in the present process in case the resonances are in different hclicity states. 

*In the first two cases, we shift the cross sections at the first and last W bins by itl.5cr of the 
systematic errors in opposite directions and make a Hnear deformation of the cross sections at the 
other intermediate bins. In the second two cases, we shift the cross sections by itl.5(T in the same 
direction at the two end bins and in the opposite direction at the center bin. We treat the largest 
observed deviation of each fit parameter from its original value as the systematic error. 
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Table 3: Resonance parameters and other variables obtained from the fit of the 
W dependence of the cross section. The first error is statistical and the second 
systematic. 



Resonance Mass Total width S Signi- 

components (MeV/c^) (MeV) (nb GeV^) ficance 

"(A) (/^(1525)) 1518±1± 3 82 ± 2 ± 3 1.22 ± 0.11 ± 0.25 > 25c7 

(B) (1.75 GeV) 1737±5±7 151 ± 22 ± 24 0.45 ± 0.09 ± 0.10 5.5cr 

(C) (2.0 GeV) 1980±2±14 297 ±12 ±6 2.63 ± 0.11 ± 0.57 > 9a 

(D) (2.3 GeV) 2327 ± 9 ± 6 275 ± 36 ± 20 0.94 ± 0.13 ± 0.28 5.3cr 

(E) a = 2.12 ± 0.17 ± 0.54 nb^/^ 

Wq = 0.04 ± 0.02 ± 0.33 GeV, b = 7.9± 1.0 ± 2.0 

Interference (A) - (E) : -1.27 ± 0.07 ± 0.12 rad 

phases (C) - (E) : +2.57 ± 0.14 ± 0.38 rad 



The small confidence level of the x^/ndf value is attributed to the inability of 
our fit function to reproduce faithfully the true behavior of the cross section rather 
than any systematic-error effects; nevertheless, we expect that we can obtain mean- 
ingful resonance parameters by this paramctrization. Wo note that there could be 
additional systematic shifts in the parameters for the particular resonant structure 
at 2.3 GeV due to the boundary of our available W range, W < 2.4 GeV. 

The statistical significances for the resonant structures arc also given in the 
last column of Table 3, They are derived from the square root of the goodness 
of fit difference between the fits with and without the corresponding resonance 
components for the structures at 1.75 GeV and 2.3 GeV. Each of the other two 
structures dominates the cross section in its energy region, and this technique is 
not suitable to estimate the significance. We show, instead, the significances for 
the excess of the core part of the peak that appears above the averaged level of 
whole the apparent peak structure, regarding the deviation from the level as the 
corresponding lower limit. 

The component (E) gives a significant contribution only at < 1.7 GeV in the 
above fit. An alternative form for (E) with two summed amplitudes having distinct 
tail decay parameters did not yield an improved goodness of fit nor a substantial 
change in the overall shape of the (E) component and the values of the other reso- 
nance parameters. We conclude that the cross section at energies above 1.7 GeV is 
dominated by the contribution from the four resonances that appear in our fit. 

We note that the size parameter S is proportional to the product Tyy{R)B{R 
K^K~) of the resonance and depends on its spin J and helicity A along the 77 axis 
as: 



S = &-k{2J +l)FT^^{R)B{R^ K+K-), (4) 
F=(2J+1) 5] /, r[di,{z)fdz. (5) 

A=0,2 "^0 

F corresponds to the fraction of the resonance component's cross section integrated 
over I cos 0*1 < 0.6 to the total, while f\ is the fraction of the resonance production 
rate in the A helicity state (with /o + /2 = 1)- The integrands in Eq. (5) are the 
square of the factors of the Wigner _D-functions. For a J = meson, only A = is 
allowed so that F = 0.6. We discuss the extraction of T^^{R)B{R K^K~) from 
this analysis in section 7. 
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Figure 6: The measured cross section of tlic process 77 K^K~ in the polar 
angular region | cos^*| < 0.6, compared with the fit described in section 6.1. The 
error bars are statistical only. The solid curve indicates the best fit. The dotted 
curve and four dashed curves indicate the contributions from component (E) and 
from the four resonance components, respectively, without interference. The dot- 
dashed curve is the contribution of the interference term. 
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6.2 Angulair dependence 



The angular dependence of the differential cross section has an enhancement at 
large angles (near | cos^*| = 0) for the lower values of W, but has a tendency of 
a forward enhancement at the higher W values. The dependence is almost flat in 
the vicinity of 2.3 GeV. We first plot the ratio of cross sections over the restricted 
range | cos^*| < 0.3 and the measured range (| cos9*\ < 0.6): 



Cll^ld\cos9*\ 

to draw out this gross behavior of the angular distribution. The W dependence of 
this ratio is plotted in Fig. 7 for W points above 1.46 GeV, where we have complete 
angular data over the range | cos 6**1 < 0.6. It has a large value at around 1.5 GeV 
and 2.0 GeV, respectively. These positions correspond to the two larger resonant 
structures designated (A) and (C) in section 6.1. 

We make a modcl-indepcndcnt angular analysis for the angular distribution of 
the K'^K~ events. The angular distribution is fitted by the sum of the partial-wave 
contributions, independently at each W bin. For the total angular momentum (J) 
of the 77 system, wc assume that the partial waves with J > 4 arc negligible. The 
odd- J waves are forbidden by the symmetric nature of the initial state and parity 
conservation. In the 77 cm. system, the differential cross section can be written as 



^ = |W(0,0) + W(2,0)|2 + |W(2,2)|2 (7) 

where 7i{J, A) is the partial-wave amplitude for the (spin, helicity)=( J, A) state, 
and is the solid angle. 

Equation (7) includes four real parameters — the sizes of the three partial-wave 
amplitudes and one relative phase between the two A = waves — while only three 
are determined independently by this measurement. (This formula is a quadratic 
function of cos^6'*.) Thus, we choose an alternative parameterization as follows. 
With aj\ = J |H(J, X)\'^dQ, as the positive definite squared partial wave amplitude 
over the whole solid angle, we parameterize the differential cross section in terms 
of three cross section-like parameters CTo, ah, and CTc and the factors of the Wigner 
D-functions as [12] 

^^^^ =aa + 5at,{dloicose*)y + 5ac{dlo{cose*)y . (8) 

The coefficients aa, ab and CTo coincide with aoo, CT20 and 0-22, respectively, if there 
is no interference between the two A = components. When we account for this 
interference, the correspondence is modified [12]: 



CTa = c^oo + Si (9) 
CTb = 0-20 + 7S/ (10) 



o-c = a22 - -zSi (11) 


Si = \/5a-ooa-2o cosV), (12) 

where tp is the relative phase between the interfering amplitudes W(0, 0) and 7i(2, 0) 
that we cannot determine independently. We note that Ua, ab and Uc can be negative 
in the presence of this interference. 
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Figure 7: The W dependence of R, the ratio of the cross sections in the angular 
ranges | cos^*| < 0.3 and | cos^*| < 0.6. 

The fit values for aa, (Tb and (Tc in each W bin are summarized in Fig. 8. The 
corresponding fit curves are shown in Fig. 4, which indicates a good parametrization 
of the measured magnitude and polar angular dependence. 

In the fits, we constrain the best-fit values of aa, at, and CTc to give a non- 
negative differential cross section throughout the full angular range, |cos0*| < 1. 
Only the fit to the lowest-energy bin, at W = 1.40 - 1.44 GeV, is affected by this 
constraint {da/d\ cos9*\ = near | cos6'*| = 0.67). 

7 Discussion 
7.1 /2'(1525) 

The resonance parameters obtained by the present measurement are summarized 
and compared with those from previous experiments in Table 4. 

Our best fit mass and total width of f^{1525) are in agreement with those from 
previous experiments [4]. In the range 1.5 - 1.6 GeV, the parameter CTc has a 
large peak. This feature is consistent with the previous determination that the 
/2(1525) is a spin-2 meson. We confirm that its production in two-photon collisions 
is dominated by the A = 2 helicity component. 

From each size parameter S, we extract Tjy{R)B{R K^K~) for the corre- 
sponding resonance R of assigned spin and helicity using Eqs. (4) and (5). These 
results are also shown in Table 4. Their systematic errors are simply scaled from 
those in S; they do not incorporate any effects from the assumptions in resonance 
formulae, background shapes and interference eff'ects. 

When we assume a pure (J, A) = (2, 2) state for /2(1525), we obtain r~^j{f2)B{f2 
-> K+K-) = 28.2 ± 2.4 (stat.) ± 5.8 (syst.) eV. This result is slightly smaller than 
the world average [4] (where isospin invariance is assumed) but is still larger than 
the ARGUS result where the interference effect is taken into account [1]. 

We note that our systematic error does not account fully for the ambiguity in the 
interference effect with the low-lying resonances /2(1270) and a2(1320), since this 
eff'ect cannot be clarified by the present data alone. We have tried an alternative fit 
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where the mass, width and size parameters of the /2(1270) and a2(1320) are fixed 
to their accepted values [4] and they are assumed to interfere with a zero relative 
phase; the quaUty of this fit is poor and the cross section is significantly larger than 
the present measurement in the mass region of the /2(1525). 

7.2 The structure around 2.0 GeV 

We see a broad resonant structure around W = 2.0 GeV in the cross section. How- 
ever, there is no remarkable structure in the individual aa, <Tb and CTc distributions 
near 2.0 GcV, nor in the total cross section crtot = o^a + crb + (Jc (wherein the inter- 
ference term si drops out). An apparent lack of consistency between this situation 
and the W-dependence of the measured cross section arises mainly from the fact 
that the variation of the parameter (T5 dominates the feature of (Ttot with its large 
contribution at small angles and conceals effects of the ac contribution, which is 
enhanced only at large angles. 

Our fit gives negative values of a,, in the wide W region above 1.6 GcV. This 
feature might imply a non-negligible (negative) interference contribution (s/) and, 
therefore, sizable values for both uoo and a2o- Given the expected dominance of 
only the 7i(2, 2) component in tensor-meson production, we still cannot explain the 
relatively narrow enhancement around | cos 0*\ = in the angular dependence for 
W between 1.8 and 2.1 GeV. 

Thus, no conclusive result is obtained for the spin of the structure near 2.0 GcV 
from our angular analysis. However, the large R value near 2.0 GeV shown in Fig. 7 
indicates that an assignment of J = 2 is qualitatively favored for this resonance, 
since both \n{2,0)\'^ and |W(2,2)p have a peak at cos 6** = 0. 

Our best fit mass and total width of the resonant structure near 2.0 GeV match 
those tabulated for the /2(2010) meson [4]. We obtain Tj^B = 61 ± 2 ± 13 eV for 
the assumption of ( J, A) = (2, 2). This result is similar in scale to that of /2(1525) 
and supports the hypothesis that this resonance is a qq meson with a sizable ss 
component [14, 15]. 

7.3 The structures around 1.75 GeV and 2.3 GeV 

Our best fit mass and width of the resonance structure at 1.75 GcV in our data 
are compatible with the corresponding values of a resonance in ref. [3] that the 
authors associate with /j(1710) in their measurement of 77 KgKg. These 
numbers arc excluded from the latest world-average values [4] of the /o(1710) since 
ref. [3] favors a spin assignment of J = 2 while ref. [4] has assigned it J = 0. The 
measurement of the 7r"'"7r~7r° final-state production from two-photon collisions [16] 
also finds a resonance near 1.75 GeV with a mass and width close to our values (they 
are shown as a2(1750) in Table 4). Ref. [4] cites this resonance as a2(1700), with 
a world-average value that includes ref. [16] as well as results from hadron-beam 
experiments. 

We cannot determine the spin-helicity of this structure from the present analysis 
alone, since the fraction of its contribution to the total cross section is small. The 
isospin of the resonant structure also can not be determined from the present ex- 
periment only. However, the difference in resonant structures between the K~^K~ 
and K'^Kg channels, where in the latter process no enhancement near 2.0 GeV is 
seen [3], can be naively explained by distinct interference effects of two or more 
resonances having different isospins. Along this line, interference between 7 = 
and 1=1 mesons — / and a mesons, respectively - is a plausible explanation for 
the contrasting behavior between the K~^K~ and KgKg channels between 1.6 and 
2.1 GeV. 
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Table 4: Resonance parameters from the results of this study compared to rele- 
vant previous measurements. M and T are the mass and the total width of the 
resonance in units of MeV/c^ and MeV, respectively, and T-^^B is the product of 
the two-photon decay width and the branching fraction to the K'^K~ final state 
for the resonance {r^^{R)B{R K^K^)) in units of cV. For the present mea- 
surement, T-y^B is derived from the S values for the stated (J, A) hypothesis. For 
prior observations (cited in the rightmost column), isospin invariance is assumed 
for comparison with our results, and any additional notes are given at the bottom 
of the table. 



Reson- 




The present 




Other possibly 




ance 




measurement 




relevant observations 




/2(1525) 


M 


1518±1±3 


/2(1525) 


M 


1525 ±5 


[4] 




r 


82 ± 2 ± 3 




F 


76± 10 








28.2 ±2.4 ±5.8 






40 ± 4.5 








((2,2) assumed) 






15.7 ±2.5 ±3.91) 


[1] 


1.75 GeV 


M 


1737±5±7 


/j(1710)^^ 


M 


1767 ± 14 


[3] 




r 


151 ±22 ±24 




V 


187 ±60 








10.3 ±2.1 ±2.3 






24.5 ±5.5 ±6.5 






((2,2) assumed) 
76 ± 15 ± 17 


/o(1710) 


M 
V 


1713 ±6 
125 ± 10 


[4] 






((0,0) assumed) 


02(1750) 


M 
F 


1752 ±21 ±4 
150 ± 110 ±34 


[16] 








02(1700) 


M 
F 


1726 ± 26 

256 ± 40 


[4] 


2.0 GeV 


M 


1980±2± 14 


/2(2010) 


M 


2011+^^ 


[4] 




r 


297 ± 12 ±6 




V 


2021^7 








61±2± 13 












((2,2) assumed) 










2.3 GeV 


M 


2327±9±6 


/j(2220) 


M 


2231.1 ±3.5 


[4] 




F 


275 ± 36 ± 20 




F 


231? 








22 ± 3 ± 6 
((2,2) assumed) 

161 ±22 ±48 
((0,0) assumed) 


/2(2300) 


M 
V 


2297 ±28 
149 ± 41 


[4] 



The units of M, F, and F-y-yS are MeV/c^, MeV, and eV, respectively. 



^) Using a coherent background. 

^) Spin 2 is reported to be dominant in the measurement [3]. 



21 



As noted earlier, the angular distribution is rather flat near W = 2.3 GeV. 
However, it is not clear that this is due to a spin property of the resonance we find 
in this mass region, since the observed evolution of the angular dependence with 
energy may come from the effect of hard parton scatterings, which have a tendency 
of giving a forward enhancement at high energies. 

In fact, the feature of the negative ac'mW > 1.6 GeV, mentioned in the previous 
subsection, may be simply an artifact of our having neglected the partial waves of 
J > 4, for example, due to hard scatterings with partons, which cause an enhance- 
ment at forward angles (corresponding to high-J waves) at these higher energies. 
Our fit compensates for the omission of any J > 4 partial waves with a negative 
(positive) offset in cTc (ct;,). We get reasonable fits of the angular distributions for 
W > 2.2 GeV with a sum of incoherent amplitudes, |'H(QCD)|2 _^ ^-^(^2, 2)|2, where 
we assume |7i(QCD)|^ has an angular dependence proportional to sin^^ [17]. 
Then, the drastic change of the angular distribution observed at 2.1 - 2.3 GeV 
provides very important information about the perturbative-QCD nature of this 
process. 

We cannot find any evidence of a narrow-width structure corresponding to the 
glueball candidate /j(2220) that is reported in the radiative decays of J/i/j near 
2.23 GeV [5] with F < 30 MoV. Howcivcr, this narrow width is not well established, 
so it is possible that a part of the wide structure that we see in this mass region 
comes from the /j(2220). We note that the identity of /j(2220) in relation with 
other resonances in the same mass region — /2(2300) and /4(2300), for example — is 
not yet clarified. 

Our size parameter 5 for the resonant structure gives T^^B = 22±3±6 eV for the 

spin-helicity hypothesis of (J, A) = (2,2). This does not contradict the upper limit 
value for /j(2220) obtained from previous two-photon measurements — for example, 
the 95% C.L. upper limit for the ifgifg channel [18] is rT,.^(/j(2220))S(/j(2220) ^ 
KgKg) < 1.0 eV for a pure (2, 2) assumption — since those analyses assume a narrow 
width of r(/7(2220)) = 20 - 30 MeV. Our analysis prefers a width about ten times 
larger. 

Working with the hypothesis of a narrow- width /j(2220), we extract an upper 
limit for its T^^B by fixing its mass and width at 2231 MeV/c^ and 23 MeV, 
respectively [4], and assuming pure (J, A) = (2, 2) production. We fit the invariant- 
mass distribution of the signal events in | cos 9* \ < 0.6 with a sum of a second order 
polynomial plus a Breit-Wigner function in the mass range between 2.13 and 2.33 
GeV/c?. We see no significant excess from the smooth polynomial level. Our 95% 
C.L. upper limit is r.^^(/,7(2220))^(/j(2220) ^ K+K') < 0.60 eV, corresponding 
to a fit including this resonance whose exceeds that of the best fit by (1.64)^. We 
account for the systematic error in the measurement by inflating the upper limit by 
1(7 of the total systematic error. This limit is an improvement on the value from 
the K'gK^ measurement which is cited from ref. [18] in the previous paragraph, 
assuming we can compare the two by isospin invariance. 

8 Conclusion 

The production of K^K~ in two-photon collisions has been studied using a 67 fb~^ 
data sample. The 77 center of mass energy dependence of the cross section for the 
process 77 K^K~ is measured in the range 1.4 - 2.4 GeV with much higher 
statistical precision than was achieved previously. 

A clear peak for /2(1525) is observed, as reported by prior experiments. We find 
three new resonant structures in the vicinities of 1.75 GeV, 2.0 GeV and 2.3 GeV, 
respectively. The mass, width and size parameters for these resonant structures 
have been obtained from the fit of the energy dependence of the cross section. 
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The structure around 2.0 GeV has sizeable contribution to the cross section, and 
the angular distribution has a large-angle enhancement in this mass region. 

The angular dependence of the differential cross section has a drastic change at 
2.1 - 2.3 GeV; below this energy, the differential cross section is more enhanced at 
large angles, and it is more enhanced at small angles above this energy. 

We do not find any signature of a narrow (F < 30 MeV) structure in the vicinity 
of 2.23 GeV. 

We hope some questions raised by this paper will be clarified by further data 
processing as well as combined analysis of different final states. On the other hand, 
better phcnomonological models taking into account the interfering resonances to- 
gether with QCD nonresonant contributions are certainly needed. 
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